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Abstract 

In this paper we prove a priori and a posteriori error estimates for a multiscale numerical 
method for computing equihbria of multilattices under an external force. The error estimates 
are derived in a W^'°° norm in one space dimension. One of the features of our analysis is 
that we establish an equivalent way of formulating the coarse-grained problem which greatly 
simplifies derivation of the error bounds (both, a priori and a posteriori). We illustrate our 
error estimates with numerical experiments. 
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1 Introduction 

Multiscale methods for modelling and simulation of microscopic features in crystalline materials 
have been very attractive to researchers of material sciences and applied mathematics in past 
two decades. In these modelling methods it is assumed that there is an underlying atomistic 
model which is the "exact" description of a material associated with certain lattice structure. 
Direct atomistic simulations using the "exact" model may not be feasible because of its huge 
number of degree of freedoms. The quasicontinuum (QC) approximation is a popular method 
to dramatically reduce the degrees of freedom of the underlying atomistic model. It was put 
forward in [29j for a simple lattice system and in [30] for a complex lattice system. Besides 
extensive application of the QC approximation in practical material simulations, there have 
been growing interest in rigorously analyzing the convergence of the QC approximation or the 
error between the "correct" and the "approximate" solutions, see, e.g., fLOl \TT\ \T3 \ [T6 l [TB I \T9 \ 
ED [22l ESI ESI ESI EU [32] , as well as a number of works attemping to design more accurate 
coarse-grained algorithms, see, e.g., [T7[ [201 [261 [271 [28] . However, most of the works, with the 
exception of |13] and [32], are for crystalline materials with a simple lattice structure. 

In this paper we consider a problem of equilibrium of an atomistic crystalline material with 
a complex lattice structure. The essential step in reducing the degrees of freedom is to coarse- 
grain the problem. The QC is one of the most efficient methods of coarse-graining the atomistic 
statics. The idea behind the QC is to introduce a piecewise affine constraints for the atoms 
in regions with smooth deformation and use the Cauchy-Born rule to define the energy of the 
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corresponding groups of atoms. To formulate the QC method for crystals with complex lattice 
(for short, complex crystals) one must account for relative shifts of simple lattices which the 
complex lattice is comprised of [30j. Our approach to model complex lattices is the framework 
of discrete homogenization, developed in our earlier paper f2j. 

We note that the idea of applying homogenization to atomistic media has appeared in the 
literature [71 El El E] • We also note that the method considered in this paper is essentially 
equivalent to the QC for complex crystals, being put in the framework of numerical homoge- 
nization [3]. However, the rigorous discrete homogenization procedure and related numerical 
method allow us to derive error estimates for the homogenized QC method, when compared to 
the solution of discretely homogenized atomistic equations. It also allows, by a reconstruction 
procedure, to approximate the original full atomistic solution. To the best of our knowledge, 
such error estimates are new. As in many numerical homogenization techniques for PDEs, there 
is no need for our numerical approximation to derive homogenized potential before-hand, since 
the effective potential is computed on the fly (see, e.g., [UlIS]). In addition, we note that the 
jji gj-j-Qp estimates in our earlier unpublished paper [2] are derived in one dimension for linear 
nearest neighbour interactions. In this paper we consider fully nonlinear multi- neighbour inter- 
actions which are technically much more difficult. Further, we will derive 14/^^'°° error estimates, 
which are more suitable for nonlinear interaction and are technically harder than those in the 

norm, and are rarely obtained even in the simple lattice case (the only estimates in W^''^ 
norm that we know of are |18| I24j). Also, we remark that we establish an equivalence of the 



coarse-grained homogenized model and the atomistic homogenized model (Lemma 4.1), which 
significantly simplifies the error analysis of the QC method. Finally, we derive both a 

priori and a posteriori error estimates. 

The regularity results of this paper are similar to those in [T3]. The main difference is 
that our results do not require a very high regularity of the external forces that was assumed 
in |13j (where, essentially, the highly smooth external forces were necessary for using inverse 
inequalities to get a W^''^ convergence from an //^-stability). 

Another related homogenization approach is the F convergence (see, e.g., [HE]) which is an 
excellent technique of finding the effective macroscopic energy from the microscopic interaction 
law, but does not yield the rates of convergence of the minimizers of the microscopic model and 
the homogenized model. 

The paper is organized as follows. In Section [2] we formulate the multiscale method for 
multilattices and state our main assumptions. In Section [3] we prove the inf-sup condition and 
regularity for the atomistic and the homogenized equations. In Section |4] we prove convergence 
of the approximate solutions to the exact ones. Finally, in Section [5] we present numerical results 
that support our analysis. 

2 Method Formulation and Main Results 

In this section after introducing the principal notations used throughout the paper, we recall 
the equations for the equilibria of multilattices and describe our multiscale numerical method. 
We then state our main convergence results. 

2.1 Atomistic Displacement and Function Spaces 

We consider an (undeformed) lattice of N atoms, £ = {e, 2e, . . . , Ne}, repeated periodically to 
occupy the entire eZ. The positions of an atom x £ C in the deformed configuration is x + u{x), 
where u = u{x) is the displacement. We will consider only eiV-periodic displacements, i.e., such 
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that u{x + eN) = u{x), thus effectively reducing the system to a finite number of degrees of 
freedom. For convenience we choose e= j^. The space of eA'^-periodic functions is denoted as 

UiC) = {u:eZ^R: u{x) = u{x + eN) Vx e eZ}. 

and its subspace of functions with zero average as 

U#ij(:) = {ueU{jO): {u)c = 0}, 

where the discrete integration (averaging) operator (•)£ is defined for u G U{jO.) by 

{u)c ■■= J^^u{x). 

We sometimes also use the notation {u{x))xec for {u)c. Also, for u,v & U{C) we define the 
pointwise product, uv, by 

uv{x) = u{x)v{x) Vx G eZ, 

and the scalar product 

{u, v)c •■= {uv)c = u{x) v{x). 

We will only consider displacements u G U^{C) since for more general displacements u{x) = 
Fx + u{x), with F G M and u G U^{C), we can adsorb Fx into the reference positions as 
u{x) = (a; + Fx) + u{x) and rescale the spatial coordinate as x = a; + Fx. 

For u = u{x) G U{jC.) we introduce the r-step discrete derivative {r eZ,r ^ 0), 

re 

For r = 1 the forward discrete derivative D^^iu we will sometimes simply be written as D^u. 
In addition to differentiation operators, we also define for u G U{jC), the translation operator 

Txu(x) := u{x + e). 

Then the r-step translation (r G Z) can be expressed as a power of T^, T^u{x) = u{x + re). 
Finally, introduce an averaging operator, 

1 ""^ 

^-,r:=-^T^, (rGZ, r>0) 
fc=o 

so that wc can write D^^r = Ax,rDx (r > 0). 

On the function space U{jC) we define the family of norms 

ll^ll"? (d'^r)/^)^^^ (1 < < oo), and ||n||oo := max |u(a;)|, 

and seminorms 

|'"|m,g := (1 < 9 < oo, m G Z, m > 0). 

The seminorms |«|m,g are extended for negative m as 

\u\m,q-= SWp{{u,v)c ■■ V E U#{jC,),\v\-rn,q' = ^ 

(1 < g < OO, (g')^ + 9^ = 1, m G Z, m < 0). 
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Note that \u\m,q are proper norms in U^{C) for all m G Z. Hence we denote spaces U{C) and 
U#{C), equipped with the respective norms, as U^''^{C) and U!^''^{C). 

We will also work with the lattice V = {1,2, . . . ,p}. For lattice functions tj = r]{y) G U(V) 
we define the operators {Dy, Dy^r, Ty, and Ay^r) and the norms similarly to functions in U{V), 
noting that the lattice spacing of P is 1 whereas the lattice spacing of £ is e. 

For functions of two variables, v = v{x, y) G fi^U(V), we will denote the full derivatives, 
translation, and averaging, by T := TxTy, Dr := ^(T^ ~ ^ •= -^i) -^r '■= \ X^fc=o^'^- Notice 
that the variables x and y are not symmetric in the definition of derivatives. If a function 
does not depend on y then the full derivatives coincide with the derivatives in x (likewise for 
translation and averaging). Hence, for functions of x only, we will often omit the subscript x in 
the operators D^, T^, A^^r- 

The following lemma, whose proof is straightforward and will be omitted, collect the useful 
facts about the above operators 

Lemma 2.1. (a) For any v £ hl[C), r G Z, r > the following estimates hold: 



\DrV\\q = \\ArDv\\q < \\Dv\\q, 
r-1 



\Dj.V — Dv\\q 



^kDkDv < ^e(r-l)p\||g. 



fe=i 1 

(b) For any v G U{jC) (8) UiJ^), r E 7., r > the following estimate holds: 

r-1 



^y,r^\\q 



k=l 



< le{r-l)\\D^v\\q. 



(2.1) 
(2.2) 



(2.3) 



2.2 Atomistic Interaction and Equilibrium 

The energy of interaction of two atoms, x e C and x + er e JC depends on three variables: 
the distance u{x + re) — u{x) between atoms x and x + er, and their positions in the reference 
configuration that are needed to account for different species of atoms. We denote such energy 
using a family of functions ^^{Dx,ruix); x) , where, for a fixed r G Z"*", is defined on (a subset 
of) M X £. The total interaction energy of the atomistic system is thus 

Ft R 
E{u) = /j2^riDru)) = /E^r(^x,rtx(x);x)\ , (2.4) 

where R is effectively the interaction radius (measured in the reference configuration). 
The equations of equilibrium are thus 

find u G s.t.: {5E{u),v)c := ^^E{u + tv)l^^ = {f,v)c eU#{C), (2.5) 

where / G U^{C) is an external force. Here 5E : U{C) — > U{C) is the Gateaux derivative of 
E : U{C) M. 



2.3 Multilattice and Homogenization 

The atoms C are assumed to be of p different species located periodically on C, and we assume 
that N G pi,. We index the atom species with V = {1,2, ... ,p}. Note that a lattice functions 
"H = ^(y) € U{V) can be related to a lattice function rj^ = r]{x/e) G U{C). 
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We define $r on an open subset of R x "P as y) := ^ri'^^v) ^ fixed r. Due to 

periodicity of the microstructure, the dependence of on y is assumed to be p-periodic, i.e. 



^riz;*) G ^i'P) for all z. For convenience of notations (e.g., in ( |2.4[ ) or (2.6)), we further 
identify, for a fixed r, the family of p scalar functions ^r{*',y), y & V, with the function 

■■ U{V) DU ^ U{V) by identifying ^r{w{y),y) with [^r{w)]{y). (Here U{V) DU ^ U{V) 
denotes a function from an open subset U of U{V) with values in U{'P).) 

We apply a homogenization to the atomistic energy to average out the microstructure; more 
precisely, to average out the dependence on y £ V. The homogenized interaction (see [2] for 
the details) is defined by 

R R 

^''{z):=J2iMz + Dy,rx{z;yy,y))yev = J2^Mz + Dy^rX{z)))r, (2.6) 

r=l r=l 

where for a fixed z G M, x(^) S hlj^{V) solves the micro problem 

R 

Y,{5Mz + Dy,,x{z)), Dy,rV)v = Vr? G U#{V), (2.7) 

r=l 

and 5<^riz;y) = £<^>riz;y). 

The homogenized interaction energy is ^'^{-^u'^)dx, whose discretized version is E^{u'^) := 
{^^{DxU^)) c- This leads to the homogenized equilibrium equations of the form 

findu° G^/# s.t.: {5<^>'^iDxu'^),Dxv)c = {f,v)c V?; G ZY#(£), (2.8) 

or, written in a strong form, 

find nO G U# s.t.: - I?,[,5$°(L',nO)] = TJ, (2.9) 

where Dx := To derive (2.9) we should use D^. = —T~^Dx- 

To extract the microstructure from the homogenized solution vP , define the corrector 

u'^ix) := T#[u^{x) + ex\DxU^{x)]x)), where (2.10) 
X#u := u- {u)c (2.11) 

and x''{z; x) := x{z', f )• Application of in the definition of u'^{x) is done for convenience so 
that u'' eU^iC). 



2.4 HQC Formulation 

Define a triangulation of the region (0, 1] by introducing the nodes of triangulation Mh C C and 
the elements Th- Each element T £ Th is defined by two nodes ^, r/ G Mh as T = £ n [^, t/), its 
interior is defined as int(r) = £ n (^, r/), and its size as hx = r] — ^. We also define the element 
size function, h G U{C), so that 

h{x) = hT Vj; G T. (2.12) 

We consider the coarse-grained spaces Uh{C) C U[C) and hlh^^{C) C U^{C) of piecewise 
affine functions. The space hlh{C) can be characterized by 

u G Uh{C) ^ Du{i - e) = Du{S) \/^eC\Afh (2.13) 

We denote the nodal basis function of V({C) associated with ^ G £ as w^, w^{x) := dx-^, 
where S is the Kronecker delta. The nodal basis function oiUh{C) associated with ^ G Mh is 
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denoted as w^{x). The functions w^, G £ \ together with w^, ^ £ J\fh, form a basis of 
IA{C). Denote the nodal interpolant Xh : U{C) — Uh{C). 

The HQC approximation to the exact atomistic problem (2.5) is 

find < G s.t.: {5E\ul),Vh)c = {F'',Vh)h yvh^Ul^iC), (2.14) 



where {•,*)h denotes the duality pairing of itii^' ^{C))* and l^h^, and G {Uf^'^{C)y is a 
numerical approximation to / G '°°(£). For convenience we extend on (Z^^' (£))* by 
requiring {F^,l)h = 0, so that {F^,I#Vh)h = {F'^,Vh)h for all Vh G (refer to ( |2.11D for the 
definition of X:^). A numerical corrector similar to ( |2.10[ ) can be introduced as follows 

ul:=I#{ul + ex'{Dul)). (2.15) 



2.5 Main results 

Before stating the main results, we introduce some additional notations. For a Banach space X 
denote Bx{xq, p) = {x £ X : \\x — xq\\ < p} — a ball centered at xq with the radius p — and call 
it the neighborhood of xq with radius p. For a mapping f : U ^ Z from an open subset U G X, 
5xf{xo) is its variational derivative at a point xq. When it causes no confusion, we may just write 
6f{xo). If / : X — 7- M with X being a Hilbert space with the scalar product (•, •)x, we identify 
6f{xo) with an element of X and write 6f{xo)x = {6f{xo),x)x', likewise the second derivative 
^'^fi^o) will be identified with a linear mapping X ^ X: {5'^f{xo)x)x' = {S'^f{xo)x,x')x- The 
space of continuous mappings f : U Z, U being bounded, will be denoted as C{U ; Z) with 
the norm ||/||c '■= ^^PxeU 11/(^)11- The space of functions whose A;-th derivative is continuous 
will be denoted as C^{U;Z) with a seminorm |/|cfc := H^'^/Hc- A space of mappings whose 
k-th derivative {k > 0) is Lipschitz continuous will be denoted as C^'^{U; Z) and the smallest 
Lipschitz constant of the k-th derivative will be denoted as | • \Qk,i. In our analysis we will 
often use the fact that if / G C'^'^"'^ then |/|cfe,i = l/lc^+i- In what follows we will express 
the statement "The quantity / is bounded by a constant that may depend on /i, . . . as 
/<Const(/i,...,/fc). 

We make the following assumptions that will allow us to apply the framework of the implicit 
function theorem (refer to Appendix [A] for its precise statement). 



Assumptions 

We assume that there exists a microstructure x* = X*iy) S U^{V) and p^ such that: 

0. The micro-deformation y + x*{y) is a strictly increasing function of y G Z. This simply 
expresses the fact that the atoms in the reference configuration are sorted by increasing 
position e{y + x*{y))- 

1. For each r G TZ and y £ V, the interaction potential ^r{*, y) is defined in a neighborhood 
U{y) C M of Dy^rX*{y) of radius p^ and $r(«,?/) G C2'i(C/(y); M). 

2. X* satisfies 

R 



J2i^'^r{Dy,rX*),Dy,rr])v =0 V?? G U#{r). 



r=l 



This assumption ensures that ex*{-) is a solution to (2.5) with / = 0. 
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3. Nearest neighbor interaction dominate: 

R 

_ Tvnn 



lmm5^^i{Dy^iX*{yy,y) - ^maic\6^^r{Dy.rX*{y);y)\ > 0. (2.16) 
y ^ y 

r=2 

Remark 2.1 (An alternative formulation of Assumption 1). It is useful to note the following 
equivalent formulation of Assumption 1 (the equivalence can he established by a straightforward 
calculation): for each r £ TZ the function '■ hl^'°°{V) DC/—)- IA^'°°{V) is defined in a 
neighborhood U of x* ^ W''^{V)) with radius p^, and G C^'^{U;U°'°°{V)). 

We next state our main results. We start with the a posteriori result. 

Theorem 2.2 (a posteriori estimate). Assume that the Assumptions 0,1,2,3 hold. For all 



G B^yi,i y{0,pf), the solution to (2.14) exists and is unique in B^i.oo{xl, Pu) ■ Moreover, 



the following a posteriori estimate holds: 

Wh ~ '"li oo < Const(cQ ^ci^'"^^) max \Du1{x) — Du^{x — e) 



+ co^||(/i-e)/||oo+ ma^ \{F\vh)t,-{f,Vh)c\. (2.17) 
Vh&Uh,#(C), 

l^/i|i,i=l 



Here C^'^' := max^gp XlT•e7^ k'^^r(», y)lco.i • 

Note that the a posteriori error estimate has a form similar to the standard FEM estimates: 
there is a term based on the jumps of the solution across boundaries of elements, a term 
consisting of summation of the external force in the interior of elements, and a term accounting 
for an approximate summation of the external force. It is worthwhile to note that for the fully 
refined mesh (i.e., where h = e), the term \\{h — e)/||oo vanishes. 

The following a priori error estimate will also be shown. 

Theorem 2.3 (a priori estimate). In addition to the Assumptions 0,1,2,3, assume that exact 
summation of the external force, i.e., that {F^,Vh)h = {f,'Vh)c- Then, for all f G B^-i,oo{0, pj). 



# 



the solution to (2.14) with the exact summation of the external force {F^,Vh)h '■= if^v] 



c 



exists and is unique in B,. i,oc{x%, Pu)- Moreover, the following a priori estimate holds: 



|'Ufe-'u|i,oo < Const(co,C^^'"^^)||/i/||oo- 

Inf-sup conditions and regularity of for the atomistic and the 
homogenized equations 



In this section we start by showing that the Assumption 3 of Section |2.5| implies the inf-sup 
conditions needed for the subsequent analysis. We then establish regularity results for the 



atomistic solution (2.5), for the micro problem (2.7), and for the homogenized solution (2 



These regularity results are essential to derive the a priori and a posteriori error estimates. 
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3.1 Inf-sup Conditions 

Lemma 3.1. Assumption 3 implies the following assertions: there exists a coercivity constant 
Co > such that the following inf-sup conditions hold 



mf sup y^(6 ^r{Dy^rX*{y);y)Dy,rri,Dy^rC) > 2co 
r?eW#(P),(^eiY#(P),^ ^ ''P 

|r;|i,cx)=l Id-, ,=1 



inf sup y^{6^<^',{Drxl)DrW,Drv)c> 2co, 

«'6W#(^).t;eW#(£),^ 

|w|l,cx)=l |i)|l,i = l 

inf sup {6^^'^{0)Dw,Dv)c = 6^^^{0)> 2co, 

h|i,oo=l |t;|i,i=l 



(3.1) 

(3.2) 
(3.3) 



where xli^) •= X*(f) '^^t? ^>*^(0) is defined by (2.6) with x(0;y) = X*{y)- 



Proof. We start with the inf-sup condition (3.2). We use the following estimate 

\{6^<^liDrXl)DrW,DrV)c\ < max\6''^',iDrXl; x)\ \\DrW\\oo \\DrV\\i 

X 

< max \6'^<l>',{Drxl;x)\ \\Dw\\oc \\Dv\\i, 

X 

for all r > 1 . For r = 1 we use Lemma 13.21 and estimate 



(3.4) 



inf sup {6^^l{Drxl)Dw,Dv)c 

k|i,oo=l |d|i,i = 1 



> i inf \\6'<!>l{Drxl)Dw 

|«'|l,oo = l 

> lmm\6^^l{DrXl;x)\. 



(3.5) 



Thus, notice that (3.2) follows from (3.4), (3.5), the assumption (2.16), and the definition 
$^(.;x) = $,(.;!). 

Proving condition (3.1) is in all ways similar to proving (3.2), with an obvious change of 
spaces U#{C) to U#{V). 



Finally, notice that (3.3) follows directly from estimating 

R 



S^q>0{0) = Y^{5^<!>r{DrX*{y);y))yer > 2co 



r=l 



using (2.16). 



□ 



Remark 3.1. The condition (3.1 ) is the same as requiring that the Hessian o/^^^ S'^^riDy,rX*] 
is positive definite, due to equivalence of the norms on finite- dimensional spaces. 



The following Lemma has been used in the proof above. 
Lemma 3.2. For u G U^, 

sup (u, Dv)c > ^ll^tlloo- 

t,ew#(£), 
l''li,i=i 



(3.6) 
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Proof. Let xi := argmax|u|. We will assume that u{xi) > without loss of generality (since 



both parts of (3.6) are invariant w.r.t. changing u to —u). Choose X2 such that u{x2) < (such 
X2 always exists for a function with zero mean) and define so that 



f 1 

2 



Dv^(x) 



X 



'2 



Xi 
X2 

otherwise. 



We obviously have |f* 



1.1 



1 and 



{U,DV^)C = \u{xi) - \u{x2) > lu{xi) = 5||U||oo- 



□ 



In the rest of the paper we will use (3.1)-(3.3) instead of using Assumption 3 directly, 



Therefore, the regularity and convergence results of this paper would hold if the ZY^'°° stability 



result (3.1)-(3.3) is proved using assumptions other than Assumption 3. Note, however, that 



the Assumption 3 is rather standard in the case of simple lattices (i.e., no dependence on y) 
and in the presence of only nearest neighbor interaction it can also be shown to be sharp. 



3.2 Regularity results 

In this section we prove our main regularity results for the atomistic and homogenized solutions. 
Instrumental for these results is a version of the Implicit Function Theorem (IFT) that we 



summarize in the Appendix (see Theorem A.l) for the convenience of the readers. For future 
use, we define 



Ccl> 



max V |^>r(»,?/)|co,i 



^e7^ 



a 



(1) 



C^^ := max max |(5^$r(», y)|co,i , and recall 



a 



(1,1) 



max > \r6^r(* 



ICO.i- 



Regularity of the Micro-problem 

Theorem 3.3. There exist pz > and p^> Q such that 



(a) For all \z\ < pz, X = x(^) satisfying (2.7) exists in C^'^(^{—pz, Pz)',U) , is unique within 
the ball U = {\\x{z) — X*|i,oo < Pxi' '^'^'^ 



IxIcM < Const(co'4')) 



(3.7) 
(3.8) 



(b) The homogenized energy density = ^^{z) is well-defined by (2.6), G C'^'^ (^{—pz, Pz)) , 
and 



(3.9) 
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Co < inf inf sup {6'^^^{z)DrW, Drv)c. 

\z\<pz |«'|l,oo=l \v\l,l = l 



(3.10) 

(3.11) 
(3.12) 



Proof. Proof of (a) 

We will apply the IFT to the mapping 



Note that (3.1) is exactly condition (ii) of the IFT. Thus, to apply the IFT, we only need to 
establish that F G C^'^ 

(2) 

Indeed, the following shows that \6^F\qo,i 
\6^Fiz',x')-6^Fiz",x")\-i,oo 

= sup I ^ Dy^^riS^-^riz' + Dy^rX) " ^^M^" + Dy,rX")]Dy,rV 



-l,oo 



< sup II ^ [6'^<^r{z' + Dy^rX) - S%r{z" + Dy,rX")]Dy,rV 
\v\l,oo=l j.^Ti 

< II ^ 5^^r{z' + Dy,rX') - S^<i>r{z" + Dy,rX") 

= max I V 5'^Mz' + Dy^rx'{v)\y) - ^^^r{z" + Dy^rXiv); v) 

< max V |52^,(.,y)|co,i|(/ - z") + Dy^rix! - x")| 



< (max \5^<^ri; y)lco.i) {\z' - z"\ + \x' - x"\i,oo), 



where we used (2.1) (and its consequence |L>rti|-i,oo < ll^lloo Vu G U{C)). The bound on 
|<5z-^lco i is obtained in the same manner. 



We hence get existence, uniqueness, and (3.8). Finally, (3.7) is obtained from |-F|co,i < C|, 
which can be proved by calculations similar to the above. 
Proof of (h ) Compute the first derivative: 



(1) 



R 



6^\z) = Yi^Mz + Dy,.x{z)), 1 + Dy,rSx{z))v = Y^^Mz + DyM^)), 



r=l 



r=l 



the last step being due to (2.7). From here we get (3.9) by taking maximum over z and recalling 
that with the assumed regularity of we have that |<I>'^|qo,i = ||5$'^||c. 
The second derivative is 



6^<^>^iz) = Yi^'^Mz + Dy^rXiz)), 1 + Dy^rSxiz))v. 



r=l 



By taking C- and 0°'^ -norms of this expression we get (|3.10l) and (|3.1ip, respectively. 

The coercivity in a neighborhood of z = 0, (3.12), is a consequence of (3.3) and continuity 

Of52«>0(z). □ 
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Regularity of the atomistic and the homogenized problems 

Define x^i^'^x) := xi^]^/^) aiid xli^) ■= X*{^/^)- We fix pz and as given by the Theorem 



3.3 



and moreover assume that py. is chosen such that < Const (1). 



Theorem 3.4. There exist pf > and p„ > such that: 



(a) For all f £ B^-i,ac{0, pj), the solution u of (2.5) exists and is unique in B^i,,x{xl, Pu)- 
# # 

Moreover, u = u{f) G C^'^{B^-i,oo{0, pf); B^i,oc{xl, Pu)), 

# # 

lli^/^illc < Cq and 
\5fu\Q0,i < Const (co, C^^^) . 



(b) For all f G By-i,oo{0, pj), the solution of (2.8) exists and is unique in Byi,oo{0, pu 

Moreover, u° = u°{f) G C^'^(5^-i,oo (0, p/); S^i,oo(0, and 

# # 



|(5/u||c < c, 



1 

' 



\6fU^\Q0,i < Const (co, C^^^) , 



and 



W'{f)\ 



2,oo 



< C, 



-1 



V/G i?^-i,co(0,p/). 



(3.13) 



In addition, the corrected solution 
(c) The following estimates hold: 



T#{u^ + ex^{Du^)) is within B^i,oc{x 



<Pu) 



,0 



""-11,00 < eConst(cQ ^C^^'"^^) |n°|2,oo- 



(3.14) 
(3.15) 



\u — u II < eConst(co^cJ,^'^^) | n°| 2,00 + e Const (p). 

Proof. Proof of (a) consists in a direct application of the IFT to (/, u) 1— t- 5E{u) — f . Assumption 
(3.2 ) guarantees the condition (ii) of the IFT; and by doing a straightforward calculation, similar 
to those in part (a) of Theorem 3.3, one can show the necessary regularity of this map. Finally, 
one should notice that (0,x*) ^ 0. 

Proof of (b). It is a standard result (cf., e.g., [23]). The proof of all the statements except 
(3.13) again consists in a direct application of the IFT to {f,u^) 1— s- 6E^(u^) — f and in all way 
similar to the proof of (a). 

To prove (3.13), we use coercivity of the homogenized problem, (3.12). For a fixed x £ C 
choose 61 G conv{Du^{x),Du^{x+e)} such that <5^>° (L>'u° (x+e) ) -5$° {Du^ (x) ) = 6'^<i>°{e){Du°{x+ 
e) — Du^{x)). By construction pu < Pz, hence S'^^^^iO) > cq > 0, therefore 



co|L>^u°(x)| < D5^''{Du''{x);x) = -Tf{x), 



0/ 



where we used (2.9), which upon taking maximum over x immediately yields (3.13). 

The possibility of choosing pu such that u'^ G B-^i,ao {xl, Pu) follows from Ixi^) — X*|i,oo < Px 



for all 1^1 < Pz which is guaranteed by Theorem 3.3 



Proof of (3.14). 



The first estimate in (3.14) is the Poincare inequality (see, e.g., |24|, Appendix A]), so we 
only need to prove the second estimate. We start with using coercivity of 6E and the fact that 
u and are solutions to (2.5) and (2.8): 



co|n^ - ^x|i,oo < \6^E{e){u'' - u)|_i,oo = \SE{u'') - 5E{u) 



-l,oo 



\5E{u^ 



5E^{u^) 



-l.OO ) 
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with some 6 G conv{n^,n} C Bj,i.oa(xl,Pu)- Thus we reduced the problem to estimating the 

consistency error, \6E{u^) — 6E^{u^)\-i^oo- 
Compute 5E{u^): 

R 

{6E{u^),v)c= Y.{S<!>l{Dru'),Drv)c 

r=l 
R 

r=l 
R 

= ^ (^T<5c|>,(D,nO(x) + Dy,rx{Du\xy, y) 

+ eD,,rTy,rX{Du\x)- y) ; y) |^^^/^, Dv{x)) _ ^ (3.16) 



r=l 



and 



V5$,(DnO(x) + Z^^,,x(^^^°(x);y); y)) ,Dv{x)) . (3.17) 

r=l 



Notice that x(-z) satisfies the equation [ ^^^^ Alj,^6(^r{z+Dy^j.x{z))\ = 0, hence X^^i ,,5$^ 
Dy^rXi^)) is constant w.r.t. y, hence 



^ /$,,(z + Dy^rxiz; y); y) 



r=l 



y=x/e 



Y,A;J^r{z + Dy,.x{z)) 
r=l 
R 

^5^r{z + Dy.rX{z)) 



V 



r=l 



Thus, combining (3.16), (3.17), and (3.18) yields 



{m{u')-bE\u''),v)c 
R 

^ [Al6<^r{Dru\x) + Dy,rx{Du\x)) + eD,,rTy,rx{Du'' {x))) 



r=l 



R 



4,,(5$.(l)n0(x) + Dy,rx{Du\x);y) ; y)]^^^,/^, I)7;(x) 



r=l 



(3.18) 



and hence \5E{u^) - 5^0(uO)|_i,oo = || T.r=i^\oo- 

In what follows we omit the arguments of DvP = Du^{x), xi^^'^^) = x{Du^{x);y), and 
(5$r(*) = 5$r(«; 2/), and likewise we omit assigning y = f before taking the ZY'^'°°(£)-norm. 

We thus estimate: 

\\Sr\\oo = \\Al5^r{DrV^ + Dy.rXiDu'^) + €D,^rTy,rX{Du'^)) 
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=: ||£:«|U + ||^:(')||oo, 



The first term is estimated as: 



= \\Al6^r{DrU^ + Dy^rXiDu'^) + eD^,rTy,rx{Du'^)) 

- Al6^r{Du^ + Dy^rXiDu"^))]]^ 

< \\5^r{DrU^ + Dy,rx{Du^) + eD^,rTy,rx{Du^)) 

- 6<!>r{Du'^ + Dy,rX{Du'>))\\^ 

< 1 5$, I co.i 1 1 15.nO - + eD.^rTy^rX {Du'^) \ \ ^ 

< |(5$r|co.i ihi"^ - l)k°l2,oo + e|xlco.i|^i°l2,oo) 

< er|5$,|co,i Const(co'4'^)l^°l2,oo, 



where in the second last step we used (2.2) and ||T, 



To estimate the term with ifi^^ we use (|2.3 



y,r ||oo 



< 1. 



< leir - l)|<5^',|co,i (|n0|2,oo + \\D,Dy^rXiDu'')\\J 

< ^e(r-l)|(5^>.|co,i (|n°|2,oo + 2|xlc".ik°l2,oo) 



< 



er|(5$,.|co,i Const (co^cW)|.,o 



2,( 



< Const (p), since 
has been estimated in ( 3.14| ). We have 



u 



Summing the estimates for £r^^ and Sr'^'^ will yield the stated result; it only remains to notice 
that C^^ < C^'^^ which implies that C^^ can be absorbed into C^'^^ 

Proof of (3.15) reduces to showing \\x''{Du^ 
ex'^{Du^) and \\u^ — u\\oo 

||x'(2;)l|oo = ||x(^)l|oo < llx(^) - X*l|oo + ||x*l|oo, 

where the first term can be estimated with the help of the Poincare inequality and Theorem 

llx - x*lloo < ilx - x*li,oo < f Px- 

To estimate the second term, recall that due to Assum ption 0, y+x*{y) is strictly increasing, 

loo < The 



3.3 



hence Dx*{y) > —1 for all y G Z, hence using Lemma 3.5 we estimate \\x* 
estimate (3.15) is thus proved. 



□ 



Lemma 3.5. Let w G V(^{V) be such that Dw{y) > —1 for all y G Then \\w\\oo 
Proof. We use the following representation of w: 

p 



< £=1 

— 2 ■ 



w{y) = ^{c-!^)Dw{y-k), 



k=l 
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which is vahd for all c G M. Choose c = 1 and estimate 



k=l 



Likewise choose c = ^ and obtain the upper bound w{y) < 



p-i 

2 ■ 



□ 



4 Proof of the main results 

In this section we prove the a posteriori and a priori error estimates. 
4.1 A Posteriori Analysis 

In order to apply our regularity results to the coarse-grained equations, we will make use of the 
following conjugate operator X^:U &s 

mw,v)c ■■= {w,Ihv)h yv,weU. (4.1) 

Note that X^tD is supported on the nodes of the triangulation Mh for all w gU, and the action 
of I"^ on w £ U can be described as distributing values of w from the interior of the intervals 
T £ Th to their endpoints. 



Lemma 4.1 (The formulation equivalent to coarse-graining). The coarse-grained problem (2.14) 
is equivalent to the following (fully atomistic) problem 

find ueU^ s.t: {6E^{u),v)c = {TlF^,v)c G U]j^{C). (4.2) 

Proof. Using the fact that the functions for ^ G i2 \ A/h, together with for ^ G AZ/i, form 
a basis oiU[C), rewrite (2.14) and (4.2) as, respectively. 



find u £U s.t. 



and 



u eUh 

{5E\u),w\)c = {f\w\)h y^eMh 
{u)c = 0, 



find n gZ^ s.t.: {5E'^{u),w^)c = {'IhF ,w^)c VCg/:\M 

{u)c = 0. 



(4.3a) 
(4.3b) 
(4.3c) 

(4.4a) 
(4.4b) 
(4.4c) 



The equations (4.3c) and (4.4c) are identical. The equations (4.3b) and (4.4b) are also equivalent 
since {I^F'^, w^)c = {F^,IhW^)h = {F'^,w^)h. It thus remains to prove equivalence of (4.3a) 
and (|4.4a|). 



Fix ^ G £ \ Mh. The right-hand side of ( 4.4a[ ) is zero, since IhW^ = and hence 

{UF\w^)c = {F\lhW^)h = 0. 



Evaluate the left-hand side of (4.4a): 



= {6E'>{u),Wi:)c = {5^\Du),Dwi)c. 
which in coordinate notation reads 

(5$° {Du{i - e)) = <5$° {Du{i)) . (4.5) 

Since is convex (cf. (3.3)), (|4.5[) is equivalent to Du{S, — e) = Du{£). Since ^ G £ \ Nh was 



arbitrary, it is further equivalent to u £Uh (cf. (2.13)). 



□ 
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Lemma 4.1 motivates us to introduce the following auxiliary problem 
find -u^"^ G U# s.t.: 



(4.6) 



We can then apply Theorem 3.4 to (2.14) and (4.6) and immediately obtain the following 
intermediate result: 

Proposition 4.2. For allZ'^F^ G i3^-i,oo (0, pj), the solution u^"^^ to (4.6) and the solution 

# 

to (2.14) both exist and are unique in B^i,oo{xl: Pu) o,nd B^i,oo(0, p^), respectively. Moreover, 

^aux(j*^/.) e C^'i and 4 = uliZlF^) E C^'i, and the 



the respective Lipschitz bounds u^"^ 
estimates 



l^aux _ 



/I- 



l.OO ) 



\Uh\2,oo < Cg ^IIX^F'' 



\ul - n^"" 



i,oo < eConst(coi4'''^)l<' 



2,oo ) 



(4.7) 
(4.8) 
(4.9) 



hold where uj^ is defined in (2.15). 



It remains to further estimate the respective quantities in Proposition 4.2 
First, we notice that e|n^|2,oo is nothing but the standard error indicator with jumps over 
elements. Indeed, for an arbitrary Uh G ^h,#, we have 



\uh\2,oo = max \D Uh{x)\ = max \D^Uh{x - e) 

xeC x&Afh 

Second, we split 



mayi \Duh{x)-Duh{x-e)\. (4.10) 



\UF^ - /|_i,oo < \iiF^'-rj''\-i,o. + - /l-i,oo 

= inax \{F'',Vh)h-{f,Vh)c\ + \nf''-f\-i,oo- 
\^h\i,i=i 



(4.11) 



Here the first term indicates how well F^ approximates the action of exact force / on the finite 
element space tlh,#- We estimate the second term using Lemma 4.3 



{TU^v)c - {f,v)c = {f,Thv)c - {f,v)c = {f,lhv - v)c < \\{h - e)f\Uv\i, 
Lemma 4.3. 



(4.12) 



where h = h{x) is defined by (2.12). 
Proof. We have 

{f,v-Ihv)c= f{x)[v -Ihv]{x) 

- e 2Z |[^ "^^"l^^^l- 

Fix T G 7/i, let ^ and rj < rj) be the two endpoints of T, and estimate, for ^ < x < r], 

\[v - Ihv]ix)\ = \vix) - ^vi^,) - ^v{'n)\ 

l{v{x)-v{0)-'^^{v{v)-v{x))\ 
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< \v{x) — v{^)\ + \v{ri) — v{x)\ 



If X = ^ then obviously [v — Ihv]{x) = 0. 
Thus, 



x'GT 



{f,v-Ihv)c < e Y max|/(x)| ^ e ^ \Dvix') 
TeTh ^ xeint(r) x'eT 

= 6 Y max|/(x)|(/iT-e) 

< \\ih-e)f\\^eY El^^(^')l- 

T&Th x'GT 



□ 



Proof of Theorem 2.2. Usmg (4.9) and (4.7) we can esthnate 



^ 1,' 



< Co^|X;:F^ - /|_i,oo + eConst(coic(''')^i-° 



The proof is then completed using relations (4.10), (4.11), ( |4.12 ). 



□ 



4.2 A Priori Estimate 

Recall that for the a priori error estimate we assume the exact summation of the external force, 
i.e., that {F'^,Vh)h = {f,Vh)c- The a priori error estimate can essentially be obtained from the 



a posteriori estimate (|2.17l) using (|4.10l) and (|4.8|). We only need to estimate \ITf^\-i,oo and 



lloo (^^^ former is needed to quantify the condition X^/ G i?^-i,oo (0, pj)) in terms of /. 

# 

This is done in the following lemma. 



Lemma 4.4. 



\nf''Uoo< i/i-i, 



||x;:/'^||oo < \\\hf\U 

Proof. To prove the first estimate, we need to prove the U^'^ stability oilh- 

\Ihv\i,i < \v\i^i. 

To prove it, start with expressing 

TeTh xeT 

Then fix T G Th, let ^ and t] {S, < "H) be the two endpoints of T, and estimate 



(4.13) 



Y\DIkvix)\ = Y 



xeT 



xeT 



\viv)-viO\ 



YDv{x) < Y \Dv{x)\. 

xGT xeT 



Hence (4.13) follows. 
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Now we can easily estimate |X^/'*|_i^oo: 

{I*j\v)c = {f,Thv)c < \f\-i,oo\Thv\i,i < |/|-i,oob|i,i, 

hence |X;;/'^|-i,oo < |/| — l,oo- 

To derive the second estimate, we test 2^/^/^ with an arbitrary v €h(: 

{ZU\v)c = {f,Ihv)c = e < e E ^ax 1/(^)1 E I^^^^K^ 

Fix r G 7/i, let and rj < rj) be the two endpoints of T, and estimate 

e \Ihv\{x) < (2^1 + < hTilHOl + hHv)\)- 



Thus, 



{rj\v)c< Yl max|/(x)|/iTQ|/(e)l + ^b(^)l) 



< 



Ten 



WhfWoo Y \v{x)\<whfWooY\<^)\ = -e\\^f\u\vWi- 

xeNh 



xeC 



□ 



The first estimate of the above lemma means that / £ B^^i,oo{0, pj) implies If^f^ G 



S^-i,^(0,p/). 



Proof of Theorem 2.3. Follows from (2.17) using (4.10), (4.8), and Lemma 4.4 



□ 



5 Numerical Examples 

We solve numerically several model problems to illustrate the performance of HQC. We consider 



a nonlinear one-dimensional model problem (Section 5.1), followed by a two-dimensional linear 
problem (Section |5.2[). 



The aim of the numerical experiments is twofold. First, we verify numerically the sharpness 
of the obtained error for the ID case. Second, we confirm that the HQC convergence result 
obtained for ID is valid in higher dimensions. 

5.1 ID 



In the first numerical example we solve the problem (2.5) with the period of spatial oscillation 



p = 2 and number of interacting neighbors i? = 3. The interaction potential is chosen as the 
Lennard-Jones potential 



with the varying equilibrium distance 

In — 



1 y is even 

9/8 y is odd. 
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Figure 1: Strain Du{x) of the solution of the ID Hnear problem: the schematically shown 
complete solution (left) and the closeup of the micro-structure for 31 atoms (right). 
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Figure 2: Results for the ID problem: error of the post-processed HQC solution u^. The error 
behaves in accordance with Theorem 12.31 



The number of atoms is = 2^^ = 16384, and the external force is taken as 

f[x) = 50sin(l-F27ra;). 
The (microscopic) strain Du{x) for such problem is shown in Fig. [Tj 

Figure [2] is aimed to illustrate that the estimate in Theorem |2.3| is sharp. Indeed, it can be 
seen that the corrected homogenized HQC solution converges to the exact solution with the 
first order in h. 

5.2 2D 

To illustrate the 2D discrete homogenization, we apply it to the following model problem. The 
atomistic lattice is £ = (0, 1]^ n e7? with e = 1/A^, the atomistic energy is 

where the set of neighbors is defined by 7^ = {(1,0),(1,1),(0,1),(— 1,1)} (we omit the neighbors 
that can be obtained by reflection around (0,0)) and the interaction coefficients as 



Such material is illustrated in Fig. [3| 



ki yi + y2 is even 
k2 yi + 2/2 is odd. 
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Figure 3: Illustration of a 2D model problem with heterogeneous interaction. 
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Figure 4: Atomic equilibrium configuration for Ni = N2 = 64 for the 2D test case. Deformation 
of the whole material (left) and a close-up (right). 



This example was motivated by the study of Friesecke and Theil [l5] , where a similar model 
was considered. Friesecke and Theil considered the model with springs similar to the one 
illustrated in Fig. [Sj which however was nonlinear due to nonzero equilibrium distances of 
the springs (so that the energy of the spring between masses xi and X2 is proportional to 

— 3;2p — Iq, where Iq is the equilibrium distance). They found that with certain values of 
parameters the lattice looses stability to non-Cauchy-Born disturbances and the lattice period 
doubles (thus the lattice ceases to be a Bravais lattice). 

The results, given with no details of actual derivation, are the following: The period of spatial 
oscillations in this case is (2,2). The function x has the form x = xO^j) = (— l)-^^^-'^ 4{k^+k2) ^ 
(here / is the 2x2 identity matrix). 

We set the values of parameters e = 2~^^, Ni = N2 = 2^^, ki = 1, /c2 = 2, A;3 = 0.25, and 
the external force 

f(^\ — in^,-cos(7r3::i)2-cos{7rX2)2 / sin(27rXl) \ _ I 

> ~ V sin(27rx2) ) 

where / is determined so that the average of /(x) is zero. The total number of degrees of 
freedom of such system is approximately 8 • 10^. The solution for such test case is shown in fig. 
|4] (the illustration is for Ni = N2 = 64). 

The atomistic domain is triangulated using nodes and K = 2t^ triangles (t = 2, 4, . . . , 2^^^). 
In each triangle a sampling domain X/j is chosen, each sampling domain contains four atoms 
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Figure 5: Illustration of a 2D triangulation. 
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Figure 6: Results for the 2D test case: error depending on the mesh size h. The error behaves 



in accordance with the ID analysis (Theorem 2.3). 



(see illustration in fig. [5|. The number of degrees of freedom of the discretized problem is 2t^ . 

The error of the solution for different mesh size h (h = 0.5, 0.25, . . . , 2^^*^) is shown in Fig. 
[6j The results are essentially the same as in ID case: the method convergences with the first 
order of mesh size in the Z//^'°°-norm. 
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A Implicit Function Theorem 

The following modification of the implicit function theorem (IFT) of Hildebrandt and Graves 
(1927), (cf. Zeidler 1986, p. 150) is used repeatedly in our analysis. 



Theorem A.l. Let X, Y , and Z be Banach spaces. Suppose that: 

(i) F ^ C^'^{U ; Z) for a neighborhood U C X x Y of {xq, i/q) and F{xo, yo) 

(ii) 6yF(xo,yo)~^ exists and is bounded. 



0. 
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Then there exist Px > and py > 0, such that 

(a) For each x G Bx{xo,px) there exists a unique solution y = y{x) G By{yQ,py) of F{x,y) = 
0. 

(b) y = y{x) is Lipschitz with the constant 

blco.i < 6o||5a;F||c < 6o|^|co.i> 

where 6o '■= 2||5yF(xo,yo)~^||- Note that \\SxF\\c < |-F|co,i due to the fact that F is 
continuously differentiable. 

(c) The derivative Sxy exists and is Lipschitz with the constant 

l^ylco.i < Const{bo\\dxF\\c,bo\5xF\co,i ,bo\6yF\QO,i) . 

Proof of estimates in (b) and (c). We assume the existence and smoothness of y{x) is proved. 

Denote 6o '■= 2||(5j/-F(a:o, 2/o))~^||- Since SyF{xo,yo) is continuous in the neighborhood of 
{xQ,yo)), we can assume that px and py are chosen small enough so that \\{6yF{x,y))-'^\\ < bo 
in Bx{xo,Px) X By{yo,py). 

Denote Fx := 5xF, Fy := 6yF, yx := 5xy. We then have Fx{x,y{x)) + Fy{x,y{x))yx{x) = 
for aU x £ Bx{xo,px), or yx{x) = -Fy{x,y{x))~^Fx{x,y{x)). 

To prove (b), estimate 

\\yx{x)\\ < \\Fy{x,y)-'\\\\Fxix,y)\\ < bo\\Fx\\c. 

To show that yx{x) is Lipschitz, fix arbitrary xi,X2 G Bx{xq, px), denote yi = y{xi) and 
2/2 = y{x2), and estimate 

\\yx{x2) - yx{xi)\\ = II - Fy{x2,y2y^Fx{x2,y2) + Fy{xi,yi)~'^ Fx{xi,yi)\\ 

< \\Fy{x2,y2y^[Fxix2,y2) - Fx{xi,yi)]\\ 

+ \\[Fy{x2,y2)~^ - Fy{xi,yi)-^]Fxixi,yi)\\ 

< 11^2/(3^2,^2)"^ II ||-^a:(a;2,y2) " ^^^(xi, yi)] || 

+ \\Fy{x2,y2)''^[Fy{xi,yi) - Fy{x2,y2)]Fy{xi,yiy^\\\\Fx{xi,yi)\\ 

< bo\Fx\co,i{\\x2 - xi\\ + \\y2 - yi\\) 

+ bl\Fy\co,i{\\x2-xi\\ + ||y2 -yi II) Iji^o. lic- 
it remains to notice that ||y2 — yi|| < ^o||-^a;||c \\x2 — xi\\ due to part (b). □ 
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